Characterization of Surface Electromagnetic Waves and Scattering on Infrared 

Metamaterial Absorbers 
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We report, for the first time, a full experimental and computational investigation of all possible 
light matter interactions on the surface of an infrared metamaterial absorber (MM A). Two channels 
of energy dissipation - diffuse scattering and generation of surface electromagnetic waves - are 
quantified in terms of their impact on specular absorption. The diffuse scattering is found to play a 
negligible roll in the absorption process, at least for wavelengths greater than the periodicity of unit 
cell. In contrast, surface electromagnetic waves are found to be generated for transverse magnetic 
(TM) polarized light at the operational wavelength of the MMA, i.e. Aq, and shorter wavelengths. 
Our computational results indicate that the highly lossy surface electromagnetic wave generated at 
Ao is responsible for the good angular dependence of absorption in TM polarization. Experimental 
results are supported by full wave three dimensional electromagnetic and eigenmode simulations. 
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Metamaterial absorbers (MMAs) [l| are a type of en- 
gineered material that obtain near-unity absorption of 
electromagnetic energy at a specified wavelength. The 
MMA operates by achieving an impedance {Z(u})) match 
to free space, while at the same time providing high loss 
- as determined by the imaginary portion of the index 
of refraction (712(0;)). A unit cell of the MMA often 
consists of two metallic layers (typically one patterned 
and one continuous) with an insulating layer lying in- 
between. The unit cell is subwavelength and periodic 
in the plane. The patterned metallic layer is an elec- 
tric ring resonator (ERR) [3|| and provides a Lorentz 
type resonant electric response (e(aj)). The ERR inter- 
acts strongly with the underlying metallic layer and a 
resonant effective magnetic response {fi{uj)) can be ob- 
tained. The first realization of such a structure was ex- 
perimentally demonstrated at microwave frequencies Q 
and designs were quickly shown in the THz infrared 
[1,0] and optical regimes [1]. The thickness of the MMA 
is small compared to the free space wavelength of light 
thus yielding large absorption coefficients. The above 
mentioned salient features of the MMA suggest they may 
be used to replace conventional absorbers and to enable 
novel applications, such as frequency-selective detectors 
0, emitters and spatial light modulators Q. 

Metamaterial absorbers operating at infrared and op- 
tical wavelengths utilize a simpler ERR geometry com- 
pared to designs demonstrated at THz and lower frequen- 
cies This permits the MMA to be fabricated more 
precisely, with less errors, thus producing higher yield 
and hence greater performance. The consequence, how- 
ever, is that the unit cell size (a) is not as sub- wavelength 
as geometries utilized at lower frequencies. Experimental 
results have shown that the resonant wavelength to unit 
cell size (Aq / a) is typically around 3 or 4 at infrared wave- 
lengths. Although IR MMAs are still able to obtain high 



absorptivity, with values of 98% shown [6] , it is unknown 
when the scattering {D(lli)) becomes non-negligible as 
the wavelength and unit cell size become close in di- 
mension. Incident energy may further generate surface 
electromagnetic wave s (S(ijj )), which may then be dissi- 
pated or re-radiated lll-ll4|. Thus in a normal specular 
characterization, both diffuse scattering and excitation 
of SEWs may appear as absorption. To date no investi- 
gation has accounted for all possible light-matter inter- 
action that may occur upon the surface of a MMA. Here 
we experimentally and computationally investigate the 
role of I?(w) and S{lu) to absorption in infrared MMAs. 
We demonstrate that scattering is small and not relevant 
in the absorption process, but that highly lossy surface 
electromagnetic waves are generated upon the MMA in 
TM polarization. 

The infrared metamaterial absorber sample was fab- 
ricated using thin-film processing combined with pho- 
tolithographic patterning. The first step was the 
electron-beam vacuum deposition of a 3 nm Cr adhesion 
layer followed by a 100 nm thick gold ground plane. The 
dielectric spacer layer of the MMA is a 195 nm thick film 
of benzocyclobutene (BCB) which was spin coated onto 
the ground plane and thermally cured under vacuum. 
Next, deep UV lithography is used to pattern the meta- 
material elements in the UV lithography resist, and a 3 
nm Cr adhesion layer followed by 75 nm thick gold layer 
are then deposited onto the nanopatterned sample. After 
performing lift-off of the Cr/Au, the remaining Cr/Au 
forms the MMA elements. A scanning electron micro- 
scope (SEM) image of a MMA unit cell is shown as the 
inset to Fig. 1(a). The patterned MM cross has the 
following dimensions: linewidth w = 500 nm, and side 
length L = 1.95 /im, with a unit cell size a — 3.2 /im. 
The fabricated samples have a relatively large surface 
area, with a total lateral size of 2.54 x 2.54 cm^, in order 
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FIG. 1: (a) Experimental diffuse scattering D{ljj) (dashed 
curve), the specular absorbance As (blue curve), and the total 
absorbance At (red curve) of the IR MM A. The inset shows 
an SEM image of the fabricated sample with dimensions: w 
— 500nm, L = 1.95 fim, and a — 3.2 ^m. (b) An expanded 
view of D{ijj) from 2-8. 5^m detailing the global minimum at 
6.13/im, as denoted by the vertical dashed line, (c) Experi- 
mental determined absorption coefScient of BCB. 
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FIG. 2: Experimental (left panels) and simulated (right pan- 
els) frequency dependent specular absorption for various inci- 
dent angles (6) in both TM (grey solid curves) and TE (black 
dash curves) polarizations. 



to enable the accurate characterization of the MMAs at 
high incident angles. 

We characterize the specular absorbance, defined as 
^5(0;) = 1 — R{iL>), vifhere R{lo) is the reflectance, shown 
as the blue curve in Fig. [TJ (The transmittance T{uj) 
was zero in all measurements due to the thick continu- 
ous gold ground plane.) The reflectance (unpolarized) 
is collected at near normal incidence (21 degrees) and is 
normalized by the reflectivity of a gold mirror. We find 
that the IR MMA obtains a prominent peak with a value 
of As= 99.8% at 6.14 fim. Outside of this main absorp- 
tive feature As is relatively low and there are small local 
maxima observed at 6.67, 6.96, 7.10, 7.99, and 9.57 /im 
due to characteristic absorptions in the dielectric spacer 
layer of BCB [15|, see Fig. [^c) for detail. 

In order to quantify the contribution of diffuse scatter- 
ing to the specular absorption shown in Fig. [Tl w e char- 
acterized the MMA using a sphere integrator lid. Il7|. A 
gold coated Lambertian scatterer (infragold®)[l^ is used 
to obtain the normalized diffuse scattered light (D) and 
normalized diffuse plus specularly reflected light (D + R). 



Radiation is incident at 13 degrees and measurements 
are performed at atmospheric conditions. In Fig. [T] (b) 
D{uj) is plotted as the dashed grey curve and is low and 
relatively featureless from 4-8.5 /im with values around 
2.2-2.3% As we proceed to wavelengths shorter than 
4 /im, D begins to increase and is 64.0% at 2.0/im. Inter- 
estingly the increase in D is fairly rapid and continuous 
for wavelength shorter than 4/im, except for some notable 
features occurring at 2.32, 2.60, 2.89 and 3.20 /im. Hav- 
ing characterized the diffuse scattering we may plot the 
total absorbance, (solid red curve of Fig. [Ha)), defined 
as At = 1 — {R{(^) + D{uj)) over the range of 2-8.5 /im 
[20I . It should be stressed that the two distinct measure- 
ments, As and At, are in good agreement over the range 
from 4-8.5 /im, and only largely deviate from each other 
where the diffuse scattering becomes significant. Figure 
[Ijb) shows a detailed view of D{uj) and it is found that 
a global minimum (^ 1.9%) in the diffuse scattering oc- 
curs at the exact position of the primary metamaterial 
absorption, denoted by the dashed vertical line. 

Having verified the negligible impact of diffuse scatter- 
ing on absorption we next investigate the possible exis- 
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tence of surface waves on MMAs and clarify their impact 
on the electromagnetic properties. We carry out a series 
of off-normal angle dependent reflectance measurements 
in both transverse electric (TE) and transverse magnetic 
(TM) polarizations. Experimental results for TE (dashed 
black curves) and TM (solid grey curves) at incident an- 
gles of 15, 30, 45, and 60 degrees are plotted as As and 
shown in the left panels of Fig. [2l For both polariza- 
tions the spectral characteristics oi As are similar for 
near normal incidence (15 degrees) - a main absorptive 
feature due to the MM A is observed near 6/zm and other 
identifiable BCB absorptive signatures are found. Pro- 
ceeding to larger incident angles we observe that, for the 
TE polarized case, the amplitude of the main peak drops 
noticeably from 90% at 15 degrees to 75% at 60 degrees. 
However, for TM polarization the peak in As remains 
high for all angles investigated, but notably two distinct 
signatures are revealed for larger incidence angles. One 
of these absorptive features - which we term here Mode A 

- is independent of incident angle and occurs at a wave- 
length of 4.83/im. The other feature - termed Mode B 

- is observed to shift to longer wavelengths as angle is 
increased (as noted by the red arrows in Fig. [2]). 

In order to gain insight into the various absorptive 
features observed in the experimental results presented 
above, we perform 3D full wave electromagnetic field 
simulations (GST Microwave Studio 2012). The MMA 
was simulated with dimensions identical to that of the 
fabricated sample. We use a Drude model for all metal- 
lic components [2l| , and a frequency independent dielec- 
tric constant of 2.06 -I- iO. 12 for BCB. Unit cell boundary 
conditions are assigned which enable us to computation- 
ally investigate the angular and polarization dependent 
absorption. Simulations provide the complex scattering 
parameters and we calculate the specular absorption as 
As = 1 - i? = 1 - The right panels of Fig. [2] show 

As for various incident angles and good agreement with 
experimental results is evident - notably both Modes A 
and B are also observed. 

We next discuss the implications of the experimental 
and computational results and focus first on Mode B. We 
note that Mode B only appears in the TM polarization 
and its peak redshifts as a function of incident angle. It 
is well-known that periodic structures can support spoof 
surface-plasmon-polaritons (SPPs) for TM polarization 



incident waves [22| . Excitation of SPPs requires the mo- 
mentum of incident light, ki, to match that of the SPP 
12, [S^, i.e. kspp — ki + kc, where kgpp and kc are 
k-vectors of SPP and the Bloch wave (reciprocal lattice 
vector), respectively. By re- writing ki in terms of free 
space light ko as ki — k^sinO, a strong angular depen- 
dence of such SPPs is evident. If Mode B is indeed a 
SPP, then its peak value should be a function of sin(0). 
We perform a linear fit to our experimental data, (not 
shown) and find that our y-intercept, kc = 27r/a, equals 
a value of 1.71x10^ (m~^). A linear fit to the curve re- 
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FIG. 3: Simulated electric field along the surface of the meta- 
material absorber at (a) 4.83 nm and (b) 6.14 /im for both 
TE, and TM polarizations, (c) \Ez\ as a function of posi- 
tion in the x-direction at 4.83/im (open triangles) and 6.14 
fim (open circles). The horizontal axis is in units of the free 
space wavelengths of Ao and Xa, where Xa denotes free space 
wavelength of Mode A. 



veals a value of a = 3.65 /im, which agrees qualitatively 
with the lattice spacing (3.2/im) of the fabricated sam- 
ples. Thus, we attribute the angular dependent feature, 
(Mode B), as a Bloch wave or spoof SPP, related to the 
periodicity of the metamaterial. 

The other absorptive feature, (Mode A), however is 
independent of incident angle, always occurring near 
4.83/im. Notably, Mode A seems to increase in strength 
as a function of increasing incident angle [24i] . Thus we 
computationally investigate its properties with light com- 
ing in at an incident angle of 90 degrees. A computational 
geometry consisting of 32 MMA unit cells is simulated 
and has a total length of 102.4/im, (see Fig. [3]). Along 
the z direction, an open boundary is assigned to the top 
of the bounding box while a perfect electric boundary is 
assigned to the bottom (the metal ground plane). By set- 
ting either a perfect electric or perfect magnetic bound- 
ary condition in the y direction, incident light along the x 
direction can be approximated as TE or TM polarization. 

The simulated transmission (not shown) shows two no- 
ticeable dips in TM polarization but, notably, is feature- 
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less for TE polarization, i.e. high transmission over the 
range investigated. Transverse magnetic polarized trans- 
mission minima occur at 6.17^m and 4.80/im, which can 
be correlated to the primary metamaterial absorption 
and Mode A in our measurements, respectively. In Fig. 
[3] (a) and (b) we plot the electric field at 6.17/im and 
4.80^m for both TM and TE polarizations. We first dis- 
cuss X—A.SOfiui, as plotted in Fig. Elja), which shows 
a TE polarized wave which does not couple to the sur- 
face. However, drastically different behavior is observed 
for the TM polarized wave. We find that the TM inci- 
dent wave is coupled to the surface of the metamaterial 
and both a decaying surface wave and a free space com- 
ponent are observed. The phase advance of the surface 
wave, compared to that of the free space wave, is different 
causing interference - a signature of light coupling with 
surface plasmons. Although quantitatively the TE and 
TM polarized waves show disparate behavior, our simula- 
tions further allow us to extract quantitative information 
by approximating both the wave number and loss of the 
SEW. For example, in Fig. [3] (a) the wavelength of the 
TM surface waves A^, may be determined by tracing the 
extrema of the electric field along the air/metamaterial 
interface. We find that the SEW yields a wavelength of 
about 4.72/im while the free space wavelength is 4.83/zm. 
Thus this mode is a type of slow waves with a k-vector 
1.02 times larger than that of free space light. In addi- 
tion, the propagation length (ispp) can be calculated by 
extracting the z component of the electric field as a func- 
tion of position (x direction), see open black circles in 
Fig. ^c) . From a decaying exponential fit to the plot we 
find that Lgpp is about 102 fj,m indicating that the SEW 
travels over an appreciable distance of 25x its free space 
wavelength of 4.83/im. Similar behavior is observed for 
TE and TM polarized waves at A = 6.13/zm, i.e. a SEW 
in TM polarization. We thus apply the same analysis in 
order to study the electric field surface profiles, see Fig. 
[3]Jb) . One noticeable difference is the TM polarized SPP 
at 6.13/im has a propagation length of 64.5 /im, which is 
only about half of that found for Mode A, shown in Fig. 
He). 

We may gain further insight into Modes A and B by 
performing cigcnmodc simulations which allow us to cal- 
culate the band structure for the k-vector parallel (A:||) to 
the MMA surface - see Fig. |4l An acoustic-like branch 
beginning from zero frequency is found with the slope 
equal to that of the free space light line. As can be 
seen, a band gap occurs before which another linearly 
dispersing mode is observed. In addition, two plasmon- 
like modes at 6 and 4.7/im can be observed. Although 
our eigenmode solver ignores losses, it should be stressed 
that these SPP-like modes occur at frequencies very close 
to the experimentally determined values of the primary 
absorption peak and mode A (red and blue dash curves 
of Fig. HI). 

We computationally and experimentally investigated 
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FIG. 4: Experimental frequency dependent absorption (left 
panel) and simulated dispersion relation in the fc|| direction 
(right panel) for the metamaterial absorber. The two horizon- 
tal dash curves indicate the peak positions of mode A (blue) 
and the primary metamaterial absorptive resonance (red). 



scattering and generation of surface electromagnetic 
waves on infrared metamaterial absorbers. Scattering 
was found to be negligible near the primary metamaterial 
absorption Aq, and only became appreciable at shorter 
wavelengths. Two surface modes were found to exist on 
metamaterial absorbers, when incident light is TM po- 
larized, one due to the periodicity of the unit cell and the 
other due to the effective optical constants of the MMA. 
The latter occurs at Aq, is highly lossy and is responsible 
for the good angular dependence of absorption in TM 
polarization. 
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